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Abstract 

After a short outline of the factorization and Grassmann picture of the one-dimensional (ID) 
Fokker-Planck (FP) equation, we consider a class of spatially-inhomogeneous solutions of the 
2D FP equation with symmetric 2D (super)potentials. We show that the spatial 
CN ■ inhomogeneities of that class of solutions can be attributed to underlying Grassmannian 

pseudo-degrees of freedom. Such an interpretation may also be applied to FP solutions in 

three and more dimensions. 
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1 Introduction 



Supersymmetric techniques are by now widely used in a rich spectrum of physical problems, 
covering as diverse fields as quantum gravity, quantum cosmology, particle physics, quantum 
field theory and statistical mechanics. In the latter area, Grassmann variables have been 
first considered by Parisi and Sourlas [jl] when they have cast into supersymmetric form a 
simple model of a scalar field coupled to a random external source. These intriguing variables 
have been most often used for Langevin equations with complex actions (Langevin equation 
formulation of quantum field theory M ) , but lately their application have been considerably 
extended. As remarked by Sourlas || the anticommuting variables are introduced more for 
combinatorial and other technical reasons and do not correspond to spin degrees of freedom 
as one might naively think. On the other hand, in the context of the Brownian diffusion 
in a ID bistable potential [Q] the simple Witten supersymmetric procedure || has led to a 
remarkable elegant way of computing the smallest nonvanishing eigenvalue Ai , which is known 
to characterize the relaxation rate towards equilibrium of a stochastic system. 

In this work, our aim is to show that a class of solutions of the 2D FP equation with 
symmetric 2D superpotential may be interpreted in terms of an underlying Grassmannian 
structure. The formalism is essentially that employed in quantum supercosmology |J], where 
however, as a rule, the cosmological potential is given and one works out the corresponding 
superpotential. The FP situation is exactly the opposite one. The formalism can be easily 
generalized to more coordinates. 

The organization of the paper is as follows. In the next section we outline Witten's scheme 
H for the ID FP equation § including also the superspace extension. In Sec. 3 we present a 
class of 2D FP solutions that can be traced to Grassmannian pseudo degrees of freedom and 
we end up with some conclusions. 

2 Factorization of the ID FP operator 

Over a decade ago, Bernstein and Brown provided a simple discussion of the correspondence 
between the ID FP equation with an arbitrary potential and Witten's supersymmetric quantum 
mechanics. The ID FP equation with constant diffusion coefficient (here normalized to unity) 
and potential drift is 



-—V(x t) - — — 

7 dt ' ypy dx 



V(x,t) (1) 



where U' = dU /dx is the drift force up to a sign, and 7 = 1/ksT in the case of the approach to 
thermal equilibrium whilst 7 = 1/vm, (y being the collision frequency and m the molecular 
mass), in chemical reaction kinetics for a system with two stable states, e.g., the trans-gauche 
izomerization process [jTj. For more algebraic symmetry, we introduced an overall - factor in 
equation (1), which can always be cast into the Schrodinger equation as follows. Any initial 
time-dependent distribution P(x, t) will relax at asymptotic times to the static solution 

Vstix) = const exp[— 7C/(x)] (2) 
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where const is a normalization constant. The evolution at intermediate times can be discussed 
conveniently by means of the celebrated ansatz 



V(x, t) = <p(x, t) exp 



(3) 



V(x, t) — ► V s t when t 



cp in imaginary time 



oo. It turns the FP evolution of V into a Schrodinger evolution for 
= —Hpp(p, where the FP Hamiltonian is a Hermitian and positive 
semidefinite operator. It is now easy to proceed with the factorization and the whole Witten 
scheme. We write H FP>1 = A^A with A = - y/jU' and A^ = + y/jU'. Thus, the 

FP superpotential is proportional to the drift force, and the superpartner Hamiltonian will be 
Hfp2 = AA^ . The two FP Hamiltonian partners are defined as usual 



-Hpp 



1,2 



ld?_ 

dx 2 



1 



1,2 



(4) 



with the potentials V\ 2 entering simple Riccati equations V^i = —jU ± H7". As it is well- 
known the great advantage of the supersymmetric procedure for the FP problem is to replace 
bistable 'bosonic' potentials with much simpler single- well 'fermionic' ones Q. The equilibrium 
distribution can always be written in terms of the FP superpotential as 



T ? st( x ) = const exp ( —2 / Wpp{x)dx 



(5) 



The time dependence of (p can be exponentiated <p(x,t) = <p(x) exp(— At) Q leading to the 
stationary Schrodinger equation 



1 d*tp(x) 
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dx 2 



+ [\ + V 1>2 (x)]<p(x) =0 . 



(6) 



In the following we shall pay particular attention to the spatial function (p(x). We would 
like to exploit some features of the formalism that may come out when considering it as a 
superfield. We shall make use of some simple rules of the Grassmannian calculus as given in 
Berezin Q. 

To introduce the superspace extension of the ID FP equation, one should write the FP 
supercharges as follows 



and 



Q = 4> 



Q = i, 



.Id . ^dU 

1 — o — 1~ l vl~^~ 
Y^y ox ox 



. 1 d _ . dU 

y/j dx dx 



(7a) 



(76) 



where tp = and ip = 9° in the Grassmann representation. The superspace ID FP Hamilto- 
nian operator will be 



H FP = {Q,Q} 



d 2 u 



dx J dx 2 



+ < 



, d d 2 U 
d9° dx 2 



(8) 
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where P = i-^-g^ is the FP momentum operator. 

Since we are interested in possible spatial inhomogeneities of "fermionic"-Grassmann origin 
we shall expand the (/>superfield in Grassmann variables as follows 

if(x,9°) = A + (x) + B (x)8° . (9) 

The conditions Qip = and Q(p = defining the 'ground-state' give A + = a + exp(— jU) and 
Bq = &oexp(7?7). In other words, A + is required to be 'square-integrable' in superspace, and 
therefore U(x) — > oo as x — * ±oo. Thus, one can see that the Bq component coming from the 
Grassmann variable can be discarded on physical reasons. However, had we considered two, 
or more Grassmann variables the conclusion would have been different, as it will be seen in 
the next section. 



3 2D FP solutions and Grassmann variables 



We now show that there exist solutions of the 2D FP equation with spatial components that 
can be naturally attributed to Grassmannian variables. Consider the following probability 

[a\ + gl{x) + gl(y)] exp(-2 7 C/) + [a 2 _] exp(+2 7 C/) (10) 



V 



where go and g\ are arbitrary functions of x and y, respectively. Suppose V is a solution of a 
2D FP equation. Then the term with the positive exponent is discarded as not physical. The 
solution given by equation (10) can be also interpreted as a solution of the 2D FP equation in 
superspace. 



In the 2D case, the supercharges read 

Q = V 



dU 
dqa 



and 



(11a) 



(116) 



The FP momentum operators are P^ 

r^=diag (1,1), i.e., ^=if v and ^ 
be written in the form 

H FP = {Q,Q} 



With {t/^, t/^}=t/^, where the metric is 



9 V , one will find the superspace FP Hamiltonian to 



if 



P»Pv + 



dU dU 



d 2 u 



dq^ dq v dq^dq u 
The 2D Grassmann representation of the <p field reads 

ip = A + + £o# + BxO 1 + A_9°9 1 



d 2 u 

dqndq u 



(12) 



(13) 



The ground-state amplitudes are determined by the conditions Qip = and Qip = with the 
Qs and <p substituted from equations (lla,b) and equation (13), respectively. Using the ansatz 
01 Bu 



^-exp(— jU), fi=0,l, these conditions lead to the following equation for / 

□2/ - 2 7 V 2 UV 2 f = 



(14) 
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where 

□2 = rf v P»Pv (15) 

and 

dU df 

V 2 [/V 2 / = V — -L (16) 

In the case of separable potential functions U(x,y) = Uq{x) + U\{y) and separable / 
functions, / = fo(x) + fi(y) one will obtain equation (14) in the form 

d 2 /o 2 dU df _ (d 2 f x 2 dU 1 dfi\ 
dx 2 dx dx \ dy 2 dy dy J 

which is in almost separable form. Denoting the /-derivatives by Fo and Fi, and the U- 
derivatives by Vq and V\ one will get 

OF 

- 2 7 V F = at (18a) 

and 

dFi 

- 2 1 V 1 F 1 = - ai , (186) 

where a\ is the separation constant. For a\ = 0, one gets easily -Fo = const exp(2jUo) and 
F\ = const ex.p(2jUi) . The solutions for B and E>i come out as follows B = 6 exp[7([/o — U\)] 
and B\ = 6iexp[— ^{Uq — U\)], where 60 and °i are constants. Thus, the probability will be 
the square of the probability amplitude in superspace Q 

V = [ a \ + bl exp(4 7 ?7o) + b\ exp(47?7i)] exp(-27?7) + a 2 _ exp(2 7 [7) (19) 

which corresponds to equation (10) above. 

When the separation constant is different from zero, one should write down the solutions 
of equations (18a,b), which are as follows 

F = 01 exp(2 7 ?7o) J ex.p(-2<yU )dx + 6 exp(2~fU ) (20a) 

and 

Fi = -01 exp(27C/i) J exp(-2 7 ?7i)dy + b x exp(27?7i) (206) 

from which the B coefficients are immediately obtained by multiplying with exp(— 7 f7). When 
a\ = we recover the previous results. 



4 Conclusions 

We have shown, on the base of a simple separable example, that 2D FP spatially-inhomogeneous 
solutions can be related to Grassmannian pseudo-degrees of freedom. Therefore one might 
think of them not only as a technical detail but also as having direct physical effects. One can 



5 



easily build higher-dimensional superfields corresponding to even more complicated spatial in- 
homogeneities of either non-equilibrium thermodynamics or chemical kinetics. In 3D, equation 
(13) is written as follows 

<p = a + + bx + x -^ vX cXq v + A_e°e 1 o 2 (21) 

where ^=0,1,2, and = ^-exp(— jU), C\ = J^exp(7?7) and the Euclidian metric is 
rf v = (1, 1, 1). equation (14) is transformed into CI3/ ± 27V3C/V3/ = 0, where the operators 
have the same meaning as in equations (15) and (16), except that the 3D metric is used. 
Equations (11) and (12) are preserved in form but the attached ^-metric is 3D. 
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